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ABSTRACT 

The Space Infrared Telescope Facility (SIRTF) will drift away from the Earth 
at ~ 0.1 AUyr -1 . Microlensing events will therefore have different characteristics 
as seen from the satellite and the Earth. From the difference, it is possible 
in principle to measure v, the transverse velocity of the lens projected onto 
the observer plane. Since v has very different values for different populations 
(disk, halo, Large Magellanic Cloud), such measurements could help identify the 
location, and hence the nature, of the lenses. I show that the method previously 
developed by Gould for measuring such satellite parallaxes fails completely in 
the case of SIRTF: it is overwhelmed by degeneracies which arise from fact that 
the Earth and satellite observations are in different band passes. I develop a new 
method which allows for observations in different band passes and yet removes 
all degeneracies. The method combines a purely ground-based measurement of 
the "parallax asymmetry" with a measurement of the delay between the time 
the event peaks at the Earth and satellite. In effect, the parallax asymmetry 
determines the component of v in the Earth-Sun direction, while the delay time 
measures the component of v in the direction of the Earth's orbit. 

Subject headings: dark matter - Galaxy: halo - gravitational lensing 
- Magellanic Clouds 



1. Introduction 

Over a dozen candidate microlensing events have been discovered toward the Large 
Magellanic Cloud (LMC) (Aubourg et al. 1993; Alcock et al. 1997a). They have typical 
Einstein crossing times t e ~ 40 days. Here, t e is related to the Einstein radius, r e by 
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where v is the transverse velocity of the lens relative to the observer-source line of sight, M 
is the mass, and d \, d\ s , and d os are the distances between the observer, lens, and source. If 
the lenses are assumed to be in the Galactic halo, they would appear to make up of order 
half the dark matter and have typical masses M ~ 0.4 M (Alcock et al. 1997a). Several 
lines of reasoning suggest that this interpretation is implausible. However, to date there are 
no plausible alternatives. 

The halo cannot be composed of M ~ 0.4 M Q hydrogen objects or they would burn and 
would easily be detected. If it were composed of white dwarfs, the white dwarfs themselves, 
their progenitors at high redshift, and the metals these produce would all be detectable in 
various ways (Fields, Freese, & Graff 1998). The most viable candidates for halo lenses 
seem to be exotic new objects such as primordial black holes - which just happen to have 
the same masses as the most common stars. In addition, if halo objects are causing events 
toward the LMC, they should also generate events of similar duration and frequency toward 
the Small Magellanic Cloud (SMC). However, the only two events discovered toward the 
SMC to date show significant evidence of being due to SMC lenses (Palanque-Delabrouille 
et al. 1998; Afonso et al. 1998; Alcock et al. 1997c, 1998; Albrow et al. 1998). 

On the other hand, if the LMC lensing events were due to lenses in the LMC bar/disk 
itself (Sahu 1994; Wu 1994), then the event rate should be much lower than observed 
(Gould 1995b), and their distribution on the sky should be more concentrated toward the 
bar (Alcock et al. 1997a). If they were due to a tidally disrupted dwarf galaxy along the 
line of sight to the LMC (Zhao 1997; Zaritsky & Lin 1998), they should be observable in 
surface brightness maps (Gould 1998a) and tracer populations like RR Lyrae stars (Alcock 
et al. 1997b) and clump giants (Bennett 1998). Hence, there are no compelling candidates 
for the lens population. 

The different possible lens populations have radically different kinematics, and one 
could therefore distinguish among them if kinematic parameters could be measured. This 
is not possible for most events because the one measured quantity, t e is a combination of 
the mass, distance, and speed (see eq. |l|). However, if a satellite were launched into solar 
orbit, the event would appear different from the satellite than from the Earth. From the 
difference, one could in effect measure the length of time it takes for the projected position 
of the lens to travel from the Earth to the satellite and also its direction of transverse 
motion. Since the distance between the Earth and satellite is known, one could thereby 
determine the two components of the transverse velocity projected onto the plane of the 
observer, 



~ d os 
v = — v. 
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Since v = |v| has values of ~ 50kms _1 for disk lenses, ~ 275 km s -1 for halo lenses, and 
^ 1000 km s -1 for LMC lenses, measurement of this quantity should distinguish well among 
components. 



2. Degeneracies 

Within the framework originally formulated by Refsdal (1966) and Gould (1994), 
the parallax-satellite measurement of v was subject to a four-fold degeneracy, including a 
two-fold degeneracy in v. These degeneracies can be understood as follows. For any given 
observer (Earth or satellite), the magnification A(u) is a function only of the projected 
lens-source separation u (Paczyhski 1986) 
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where to is the time of maximum magnification and (3 is the impact parameter. In units 
of the projected Einstein ring f e = (d os /d\ s )r e , the separation between the Earth and the 
satellite along the direction of the lens motion is simply At /t e , where 



At = t 



0,5 



t 



0,® 



(4) 



is the difference in the observed times of maximum. However, the separation in the direction 
normal to the lens motion (in units of f e ) has four possible values, ±A/3 + and ±A/3_, where 
A/3 ± = |/3 s ±Ab|. Since 



Au 



d 



Au = (At /t e ,A/3), 



(5) 



this engenders a four-fold degeneracy in v in direction and a two-fold degeneracy in 
amplitude. Here d is the Earth-satellite distance. 

It is, however, possible to break this degeneracy by taking advantage of the fact that 
the velocities of the Earth and satellite differ in the Af3 direction, leading to a difference in 
observed time scales, At e = 2-nyi~ l A f3t 2 e , where At e = t ei s — t e ,e, and where I have adopted 
the conventions that the satellite is trailing the Earth and that if v is outward from the Sun, 
then A/5 > 0. Hence, if this time difference can be measured, Af3 can be unambiguously 
determined (Gould 1995a), 



Af3 = 



1 At P 
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However, measurement of the time-scale difference introduces an additional degeneracy 
which must be dealt with very carefully. The observed flux from the star is actually a 
function of five parameters 

F{t- t , (3, t e , F , B) = F A[u(t; t , (3, Q] + B, (7) 

where F is the source flux and B is any background light that falls into the aperture 
but is not lensed. Measurement of B is highly correlated both with (3 and t e because the 
effect on the light curve of with changing all of these three parameters is even in (t — to) 
and very similar to one another. Hence, it is all but impossible to measure the small 
time-scale difference between the Earth and the satellite if each light curve is fitted with its 
own background parameter, B. When I proposed this method (Gould 1995a), I therefore 
explicitly assumed that the background light was identical for the two sources, B s = B®. 
Physically, this condition can be attained if the filters have identical (or nearly identical) 
transmission properties, and if the images are convolved to the same seeing. With this 
constraint, t e> s is still highly correlated with Bs, and t e> ® is highly correlated with £> ffi . 
However, since B s = B 9 , the difference At e = t ei s — t e ,® is we U constrained. 

The Space Infrared Telescope Facility (SIRTF) will be launched into solar orbit early in 
the next decade, and will drift away from the Earth at ~ 0.1 AUyr -1 . In this respect, it is in 
an excellent position to measure microlens parallaxes. However, the bluest band that it can 
observe is L band (~ 3.6 fim). Because of the complexity of atmospheric transmission at L, 
it is not possible to mimic the space-based detector response from the ground. Moreover, 
the ground-based background in L is so high that it is not practical to monitor typical LMC 
sources (V ^ 20) from the ground. Hence, one cannot insure Bs = B®, and so it is not 
possible for SIRTF to measure parallaxes using my original single-satellite approach (Gould 
1995a). 

3. A New Parallax Method 

To recapitulate, SIRTF has no trouble measuring one component of Au, namely Ato, 
which is unambiguously given by the difference in times of peak magnifications. However, it 
cannot measure the other component, A/3, either directly or through measurement of At e 
(via eq. §) because both parameters are degenerate with the blended light, B. The key to 
resolving this problem is to notice that equation (|6]) arises from the Earth-satellite velocity 
difference in the Aj3 direction. However, the Earth itself undergoes a velocity change in this 
direction during the course of the event by an amount Av ~ Qt e v^, where = 30kms _1 
is the speed of the Earth. This purely ground-based parallax effect leads to an asymmetry 



- 5 - 



in the light curve which can, in effect, be used to measure A/3 and hence complete the 
measurement of Au. 

Actually, such parallax asymmetries have a long history. I showed that for events of 
sufficiently long duration, it would be possible to measure a complete parallax from the 
ground (Gould 1992). One such parallax measurement has been published for a long event 
(t e ~ 110 days) seen toward the bulge (Alcock et al. 1995) and several others have been 
observed (Bennett 1997). For shorter events (fi©t e <^ 1), the acceleration of the Earth can 
be approximated as constant over the course of the event. Equation ([3|) is then replaced by 
(Gould, Miralda-Escude, & Bahcall 1994) 

u(t;t ,/Me,7) = 

where 

7 = ^ffi*e cos (f), (9) 
v 

cf) is the angle between v and the Earth-Sun separation, and where I have made use of the 
fact that the LMC is approximately at the ecliptic pole. Gould et al. (1994) noted that 
even for very short events, halo lenses (v ~ 275 km s -1 ) could be distinguished from disk 
lenses (v ~ 50kms _1 ), at least statistically. The problem is that in any individual case, if 7 
were measured to be consistent with zero, one would not know whether the lens were in the 
halo (where v is large) or in the disk (and cos0 just happened to be small). Moreover, the 
principal question about the location of the lenses is not halo vs. disk, but halo vs. LMC. 
For short events it would be extremely difficult to distinguish between halo and LMC lenses 
using this technique, even statistically. 

However, the typical events observed toward the LMC now turn out to be considerably 
longer, t e ~ 40 days. Indeed, their long time scale is a major puzzle. Hence, I recently 
proposed that LMC events be intensively monitored to search for this effect (Gould 1998b). 
I showed that one could distinguish statistically between the halo-lens and LMC-lens 
scenarios by observing 15-30 events over 5 years. The GMAN collaboration (Alcock et al. 
1997d) is routinely monitoring LMC events, but probably not intensively enough to detect 
this effect. In any event, the fact that one measures only the degenerate combination, 
v sec 0, and not the two parameters separately, means that one cannot resolve the nature of 
the events on an individual basis. 



t-tr 
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3.1. Overview 

The new approach is to combine a ground-based measurement of 7 (and of course to,©) 
with a space-based measurement of io,s (and therefore At ) to completely determine Au. 
With this is mind, I write u in a basis that is rotated relative to equation (|||), 

Au = (Au x , Au y ), Au x = — cos 9 + A/3 sin 6, Au y = -— sin6» + A(3cos9 (10) 

Here, Au x and Au y are anti-parallel and outward normal to the direction of the Earth's 
motion at the midpoint of the event, while Ato and A/3 are anti-parallel and normal to the 
satellite-Earth separation vector. The rotation angle 6 is defined by 

SiD ^2XU' < U > 

This is convenient because Au y can be simply expressed in terms of observables (see eqs. || 
and§, 

AUy = M ( n ©^©)" V (12) 

From both a conceptual and a practical point of view, it is useful to think of the 
measurement process as first determining the four Earth parameters, t ,e> /$©> te,®, and 7®, 
and then using these to predict the four analogous satellite quantities as a function of the 
(unknown) parameter Au x . 

to,s = to,© + At , j^- = Au x cos6 - 2(fi e t ei0 )~ 2 7 e sin 2 6 1 (13) 

/3 S =|/3©±A0|, A/3 = A M;E sin^ + (fi ffi te,e)" 2 7e sin 2^ (14) 

t e ,s = *e,© + At e , ^ = ^u x Q e t e ^ sin 6 1 + (fi ffi t e) ©) -1 7 e sin 26* (15) 

7s = Au,.(ft e t Cie ) 2 cos0 + 7© cos 26* (16) 



The satellite measurements are then used to determine Au x . 
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Equations ( |ilj) and (15) yield very little information about Au x because (3 and t e 
are poorly determined, being strongly correlated with B. In addition, equation ([14]) is 
ambiguous because, in the limit fi©t e / <; 1, the ground-based observations do not yield 
information about whether the lens passed on the satellite side or the opposite side of 
the Earth. Thus, equation fllil) contains essentially no information about Au x . It is not 
immediately obvious, but I show in § 3.2 that equation (|16D also has relatively information 
about Au x . Hence, in agreement with our naive expectation, Au x is mainly determined by 
measuring Ato- However, from equation ([13]) we see that uncertainties in the measurement 
of 7® can propagate into the Au x determination. 



3.2. Ground-based Observations 

I have previously discussed in some detail the problem of early identification and 
intensive monitoring of microlensing events toward the LMC (Gould 1998b). As in that 
paper, I evaluate the covariance matrix Cjj of the six parameters <2j(= to, P, t e , 7, F , B) 
specified in equations (|]), (|7|), and by considering a series of measurements at times tk, 
and with errors a^, 

,-1 , ^ - 2 9F(t k ) dF(t k ) 

After taking the derivatives dF(tk)/da>i, I evaluate them assuming 5 = 7 = 0. I assume 
that the errors are photon limited, i.e., a k = cro-F [v4(t fc )] 1//2 . (This assumption differs 
somewhat from Gould 1998b.) I assume that these intensive observations are triggered 
when the event enters the Einstein ring (itinit — 1, ^.(^init) = 1-34) and end at t — to + 1.5t e , 
and that they are carried on uniformly at a rate per day in the interval. I then find an 
uncertainty in the determination of all parameters. In particular, for the key parameters, 
to, 7, and (3. I find 

Ut0 - ^ ;R(fi), (18) 



t e (Ate/day) 1 /?' 



and 



- = a ° ^ = rm ?Lm,2 t W> ( 19 ) 



7 (Nte/d&yy/*"^" ^ (At e /day)V 
where R((3), S(/3), T(0) are shown in Figure [l]. 



- 8- 



These results shed light on three questions regarding equations (|12|), (|T3"D , and (|T6|). 
First, the ratio of the contributions of <r 7 and a to to the error in the determination of Au x 
from equation ( |T3"D is 

2(n e t e )- 2 sin 2 0-^- ~ 0.4 ( -*—) 2 ( -^-) ~ 2 -^L. (20) 
V ® e; di /t e V0.2AUy V 40 day/ 10R(p) v ; 

Thus, for typical parameters, the two sources of error are comparable. See Figure [l]. 

Next, I note that the error in Au x induced by the "7©" term in equation ([16]) is larger 
than the error induced by the "At " term in equation flT3| ) by a factor ~ S(f3)/R(f3) ~ 10. 
Hence, for typical parameters, Au x is constrained primarily by equation flT3|). This result 
has the important consequence that the SIRTF observations should be optimized to 
constrain to,s rather than 75 or any other parameters (see §3.3). 

Finally, the error in Au y is larger by a factor ~ 2(d/AU) _1 than the error induced 
by the 7© term in equation (|1"3|). By equation (|20| ) the latter is comparable to the total 
error in Au x . Hence, the (purely ground-based) error in Au y is substantially larger than 
the ground-based contribution to Au x . That is, the limit imposed by the ground-based 
observations on the overall precision of the measurement is set by the ground-based 
measurement of Au y , i.e., 7. From equations (H), (|l"2"D, and flT9f) , the error in Au y expressed 
as a fraction of Am is 

= a -2 =0 17 N- 1 ' 2 * ( te ) ^S- (21) 

Au 7sec0 ' 0.01 275 km s" 1 V 40 day/ 8 v ; 

Thus, for a robust detection of parallax for a halo event with typical parameters requires 
oq = 1% photometry, on average once per day. 



3.3. SIRTF Observations 

There are a number of constraints that affect SIRTF observations that do not 
affect ground-based observations, the most important of which is scheduling. While the 
operational plan for SIRTF is not finalized, it is expected that the IRAC camera (used for L 
band photometry) will share time equally with two other instruments, rotating for example 
one week on, two weeks off. Thus, the type of schedule envisaged for the ground-based 
observations (exposures once or several times per day) are out of the question for SIRTF. 
Instead, observations must be concentrated in a few critical periods, each less than a week 
and separated by several weeks. This appears at first to pose major difficulties, but in fact 
such constraints are perfectly compatible with an optimal observing plan. 
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To devise an optimal strategy, first recall that the flux is a function of time and six 
parameters, F(t; Oj). The derivatives of this flux with respect to the parameters, dF/dai, are 
even in (t — t ) for four parameters (a, = (3, t e , F , B) and odd for the other two (a; = t , 7). 
Hence, any set of observations that was symmetric in (t — to) would automatically decouple 
the errors in ((3,t e , F , B) from those in (to, 7)- This would be good because, from the 
discussion in § 3.1 and § 3.2, t and 7 provide essentially all of the useful information. I 
focus first on the simplest symmetric case, two observations placed at t± = t ± rt e , where 
t is a parameter. In fact, a third observation is required to establish the baseline, Fq + B. 
I will return to this point below, but for the moment I assume that the baseline is known 
with perfect precision. 

Now, if it were literally true that there was no information about (3, t e , and F — B, then 
of course a light curve fit to only two points would be completely unconstrained. However, 
from the ground-based data and equations (|T4l) and (|15|), (3s and £ e ,s are reasonably 



constrained. I will discuss this in more detail below. Here I focus on what specifically can 
be learned about t s, and 75 from this symmetric pair of observations. The (2 x 2) inverse 
covariance matrix associated with t (normalized to t e ) and 7 is (see eqs. ||, [7], and [T7]) 

64 ( 2t 2 -t 4 \ 

MV*e,7) " u 5 (u 2 + 4 )5 /2 (u 2 + 2 )<r 2 \- T * T «/2 ) ' ^ 

where u 2 = r 2 + f3 2 . This matrix is of course degenerate between to and 7. I assume for 
the moment that the information in equation flip] ) is sufficient to break this degeneracy. 
(I check this assumption below.) The error in to/t e from these two measurements is then 
a tQ ~ [6(1, l)] -1 / 2 . This error is minimized approximately at r ~ (2/3) 1 / 2 /? at which point 

Y~{W ^ {a * = ( 5 /3) 1/4 j8 1/2 *o), (at r = (2/3) 1 / 2 /?) (23) 



where a* is the approximate fractional flux error of the two observations, and where I 
have made the evaluations using the approximation A(u) ~ u^ 1 . Thus, for typical values, 
(3 ~ 0.4 and t e ~ 40 days, the two observations are separated by ~ 26 days, and the error in 
t is cr to /t e ~ 0.6cr*. 



3.4. Correlations 



To arrive at this estimate, I have argued or assumed that one can ignore the numerous 
correlations among the 12 observable parameters (six from the Earth and six from the 
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satellite). I test these assumptions by simulating a fit based on the observations as outlined 
in § 3.2 and § 3.2. In order to do so, I must choose a relative scale of errors for the Earth 
and satellite observations. The observations should be designed so that the errors in Au x 
and Au y are roughly comparable. On the assumption that [6(1, l)] -1 / 2 in equation ( P3"D 
gives a good estimate of the error in to and that this error dominates the error in Au x 
(see eq. H), I initially assume that (25/12) 1/2 (3a* = (d/AU)(fi ffi t e )- 2 <7 7 (see eq. pD- For 



defmiteness, I initially assume d = 0.2 AU and t e ^ = 40 days. 

I fit for a total of 13 parameters including 12 observables (a^, i = 1 . . .6 for to, fl, t e , 
Fq, B, and 7 as seen from the Earth, and ai, i = 7. . . 12 for same parameters as seen from 
the satellite) plus 013 = Au x which is a derived parameter. The inverse covariance matrix 
bij is given by the sum of four types of terms. First, there are terms of the form given by 
equation ( |T7D for Earth-based observations which affect b^, i,j = 1 . . .6. Second, there 
are terms of the same form for satellite-based observations which affect 6^, i, j = 7. . . 12. 
Third there are contributions to 6^- from the constraints (|I"3"D , (1~5|), and ([16]). Each of 



these can be written in the form Yl,iCti a i = 0. For example, for equation (|i~6l), ai2 = 1, 
a 13 = — (fl^t e ^)~ 2 cos 9, a 6 = — cos 2$, and = for all other i. These constraints lead to 
contributions to 6^ of the form a^aj/Q 2 where Q is an arbitrarily small number. Finally, 
equation ( |T4"D leads to a similar constraint, except that there is a discrete uncertainty in 
the sign of A/3. Hence the constraint is less certain and so has the form aiaj/(A/3) 2 . For 
defmiteness, I choose Af3 = 5a & u on the assumption that the observations have been 
structured to detect A(3 at the 5 a level. 

I then find for a pair of observations that are exactly symmetric about to,s, that for 
(3 — 0.1 . . . 0.5, the errors are higher than my naive expectations by factors / = 1.03 1.07, 
1.14, 1.24, and 1.36. If the various assumptions that I have made in the analytic derivation 
had all held exactly, these ratios would all be unity. The deviations from unity are partly 
accounted for by the fact that to derive equation (|23|), I evaluated equation ( |22"D assuming 
(tt 2 + 4) 5 / 2 (w 2 -|-2)/64 = 1, whereas it is actually slightly higher and increases with increasing 
(3. Another factor is that in equation (||), [6(2,2)/6(l, 1)] 1/2 = /3 2 /3. Thus, as (3 increases 
the role of the 6(2, 2) (7) term is relatively less well constrained by equation (|i~6|). In brief, 
the estimate (^) for the required photometric precision is basically accurate but is slightly 
too optimistic for (3 ^ 0.3. 

However, the assumption that the pair of observations is exactly symmetric about to,s 
is too idealized. Even if one had perfect freedom to schedule the observations, and even 
if one knew £ ,®> ^e,e> an d 7© exactly from the ground-based observations, according to 
equation ( |i~3|) there would still be an uncertainty Au x t e<9 cos 9 in the predicted value of 
At (and so of to,s), where Au x is an unknown quantity still to be measured. For typical 
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parameters, At ~ d/v ~ (9(1 day). In addition, at the time that the second observation 
is planned, there will be some measurement uncertainties in t ,e> ^e,®> an d 7®- This will 
lead to an additional uncertainty in the prediction of io,s> although this uncertainty will 
probably be less than 1 day. The biggest potential problem is that it may not be possible 
to schedule the second set of observations exactly when one would like. 

Regardless of the reason, if the observations are not symmetric about to,5; then the 
errors in (3, t e , and F — B will not decouple from those in to, and 7. These three quantities 
are relatively poorly determined, so the degradation of the precision could be significant 
if the asymmetry of the observations about the peak is sufficiently large. To quantify 
this effect, I imagine that the two observations take place at t± = t + St ± (2/3) 1 / 2 /5t e . 
That is, their midpoint is displaced from t by a time St, but the separation between the 
observations is the same as in the optimal case (eq. 1531) . (The effect of using a different 



separation can easily be judged from the (1,1) component eq. ^). Since the errors are a 
minimum for St = 0, we expect that they will be quadratic in St. 

Figure shows the ratio, /, of the true error in Au x to the naive error as a function 
of St for (3 = 0.2, 0.3, 0.4, and 0.5. For definiteness, I have as before assumed d = 0.2 AU 
and t e = 40 days. The error in Au x becomes seriously degraded if St is more than about 
two days. The effect is worse for low /3 because the effective time scale t e ff = /3t e is shorter, 
so the asymmetry of the observations is more severe for fixed St. This result emphasizes 
the importance of scheduling the observations to be as symmetric as possible. However, it 
also shows that the inevitable ~ 1 day errors in estimating £o,s will not seriously affect the 
precision. 

If the satellite separation is reduced to d — 0.1 AU, but t e remains at 40 days, then the 
results shown in Figure |^ remain qualitatively the same. On the other hand, if d remains at 
0.2 AU while t e is reduced to 20 days, then all the curves rise about twice as rapidly. That 
is, the figure remains approximately accurate for all relevant values of d and t e provided 
that the abscissa is labeled u (St/t e ) x (40 days)". 



3.5. Constraining the Baseline 

The calculations of the previous section were somewhat idealized in that they assumed 
that the baseline (F + B) is known exactly. For satellite observations such exaggerated 
precision would come at a very high cost. In fact, observations of the baseline need not be 
very intensive. This can be understood as follows. Let the fluxes measured on opposite sides 
of the peak be F\ and F 2 , and suppose that the measurements are nearly symmetric so that 
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AF = Fx — F 2 <C F* = (Fi + F 2 )/2. From the ground-based measurements combined with 
equations (|I3|)-(|T6|) (plus the fact that |Au| <C 1), one knows t ,5, Ps, te,s, an d Is to within 
a few percent, and hence one also knows A at the time of the observations with similar 
precision. The quantity that gives information about Ato is 5 A/ A = 5F/(AFo). That is, 
uncertainty in the estimate of Fq will degrade the precision of Ato only if its fractional error 
is of the same order or larger than the fractional error in At - Since the latter is not likely 
to be much lower than ~ 20%, only a relatively crude measurement of F is necessary. The 
near-peak measurements give AFq + B fairly precisely, and the baseline gives Fq + B. The 
error in the difference, (A — 1)F , will thus be of the same order as the error in the baseline. 
Since A is known relatively well, Fq will be determined with similar fractional accuracy as 
the baseline. Thus, the baseline measurement can be an order of magnitude or more less 
precise than the peak measurements. 

I find numerically that if the baseline exposure time is equal to the exposure time 
for each of the near-peak observations, then the precision of the determination of Au x is 
affected by 1% or less. Even if the exposure time is reduced by a factor 10, the precision 
of Au x is degraded by 10% or less. Thus, the total required observation time is well 
approximated by the time required for the two near-peak observations. 

4. Practical Considerations 

Here I consider the problems of timely event recognition, signal-to-noise ratio (S/N), 
image analysis, telescope time requirements, and backgrounds. 

4.1. Event Recognition 

The event must be recognized sufficiently early for two reasons. First, the precision of 
the ground-based measurement of 7® (and so Au y ) depends critically on when the intensive 
follow-up observations begin. Second, the characteristics of the event (t , (3, t e , and F ) 
must be sufficiently well understood from the initial ground-based follow-up observations 
to plan the satellite observations. These latter are, in their nature, target-of-opportunity 
observations and so will inevitably require some rescheduling. Moreover, given the rotation 
of instruments, there is only a 1/3 probability that the IRAC camera will be scheduled for 
use at the optimal time for the first exposure and another, independent 1/3 probability that 
it will be scheduled for use at the optimal time for the second exposure. Thus, a complicated 
decision process will be necessary to balance the requirements of these observations with 
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other aspects of the SIRTF mission, and it is important that sufficient information about 
the event be available to make a rational decision on a timely basis. 

Early recognition presents a significant challenge for the proposed observations. At 
present, the MACHO collaboration alerts on events when they surpass magnification A = 1.6 
(wi n it = 0.75) rather than A = 1.34 (u init = 1) (A. Becker 1998, private communication). 
Figure |3] compares the S/N function S(j3) (see eqs. [19] and |21|) for these two values of M ini t. 



It is clear from this figure that the S/N is severely degraded for j3 ^ 0.45 at u init = 0.75. 
A closely related problem is that the optimal time for the initial satellite observation is at 
u = (5/3) 1 / 2 /? = 1.29/3 (see eq. |23|). Thus, for (3 = 0.5, the optimal time for the first satellite 
observation is at u = 0.65 which is only a time 0.15 t e ~ 6 days after the present-day initial 
alert. 

It is therefore important, though not absolutely critical, to alert on events earlier than 
the present standard. One approach would be to initiate aggressive follow-up observations 
at a lower threshold, and to weed out the false alerts through these intensive observations. 
This would probably require a dedicated or nearly-dedicated 1 m telescope. Another 
approach would be to obtain higher S/N during the initial microlensing search observations. 
C. Stubbs (1998 private communication) and collaborators are trying to organize a search 
with a 2.5 m telescope, 1" seeing, and a 1 deg 2 camera, which would represent a factor 7 
improvement in S/N relative to the present MACHO setup. It might well then be possible 
to alert on events at lower magnification. However, even if neither of these more aggressive 
programs are implemented, one could still carry out parallax measurements by restricting 
events to those with j3 ^ 0.45. 



4.2. Image Analysis 

I assume that the observations will be analyzed using image-differencing techniques 
which have been pioneered by Tomaney & Crotts (1996) and Ansari et al. (1997) to find 
microlensing events of unresolved stars in M31. Melchior et al. (1998), Tomaney (1998), and 
Alard & Lupton (1998) have further refined these techniques for application to photometry 
of resolved (or partially resolved) sources such as those that are routinely monitored in the 
LMC, SMC, and Galactic bulge. A version of the Alard-Lupton algorithm is now used by 
the EROS collaboration to make precise light-curve measurements for events found using 
their normal (more standard) photometry (Afonso et al. 1998). 

The application of these techniques to the SIRTF IRAC camera should be relatively 
straight forward, but not completely trivial. The pixel size of this camera is l."2 while the 
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point spread function (PSF) is about 2". This is somewhat larger than the diffraction limit 
at 3.6 /im (with a 0.85 m telescope) of ~ 1". However, the PSF will still be undersampled. 
The two exposures that are to be differenced will be separated by ~ 25 days and therefore 
rotated relative to one another by ~ 25°. Hence, the source star will not fall on exactly 
the same parts of the pixels in the two exposures. Variations in sensitivity across a single 
pixel would therefore make it impossible to obtain the required photometric precision (see 
§ 4.3) by point-and-stare observations. Hence, it will be necessary to dither the telescope 
by fractions of a pixel many times in order to enlarge the effective PSF and so smooth over 
the pixel-sensitivity variations. 

The normal practice in image differencing is to form a "template image" by combining 
several images at baseline, and then to subtract this from the event images. For the present 
case the procedure is quite different. The main information comes from directly differencing 
the two images taken symmetrically about the peak of the event. The shorter baseline 
exposure is not used as a template. Rather, the difference between this image and the 
average of the two images near peak is used to extract the less precise auxiliary information 
about F . See § 3.5. 

I take note of a minor technical point. Image differencing automatically removes the 
background sources B from the analysis. This means that a standard microlensing light 
curve is fit to four parameters rather than five. That is, one fits to 

F{t; t , P, t e , F ) = F [A(t : t , (3, t e ) - 1], (24) 

instead of equation ([I]). For this reason, it is sometimes mistakenly believed that the fit 
is less affected by the degeneracies associated with background light. Recall from § 2, 
that removal of these degeneracies is the central problem addressed by this paper. Why is 
equation (|24|) no more constraining than equation (|7D? The latter can be formally rewritten 

F{t; t , P, t e , F Q , B) = F [A{t : t , (3, t e ) -1] + B, (25) 

where B = F Q + B is the baseline flux. In normal microlensing observations, B is extremely 
well determined because there are a large number of observations at baseline. Hence, it is 
effectively decoupled from the other parameters in equation (p5|). On the other hand, F in 
equation (p4"|) refers to the flux above baseline. That is, equation (|24T) implicitly assumes 
that the baseline is well observed and that an essentially perfect template has been formed 
from these observations. Thus, in reality, equation (pl|) has the same information content 
as equation (^) or equation (p5|) 
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4.3. Telescope Time 

I have already discussed at some length the requirements for the ground-based 
observations to measure 7 e and thus Au y (Gould 1998b). The basic result is summarized 
by equation ([H]). I focus here on the requirements for SIRTF observations. 

In § 3.4, I introduced the quantity a*, the fractional photometry error for 
symmetric pair of observations near the peak. Recall that the error in Au x is given by 
<Jau x — (25/12) 1//2 //3cr,,, where / is the correction factor shown in Figure 0. For reasonably 
symmetric pairs of observations, 1 < / < 1.36 The companion to equation (|21|) is then 

^'V.^^o.w^A^i^^f^)- 1 . (26) 



Am V127 Jr d ' J 0.01 0.4 275kms- 1 40day V0.2AU 

Thus, for a robust detection of parallax for a halo event with typical parameters requires 
o"* = 1% photometry for two each of two observations near the peak, plus a shorter exposure 
of the baseline (see § 3.5). This is to be compared with the ground-based requirement of 
1% photometry, once per day for a period ~ 2.5 t e . 

The IRAC detector records 0.7 e~ s _1 /iJy -1 , and the sky plus dark current is expected 
to be 3 e~ s _1 pixel -1 (J. Hora 1998 private communication). The pixel size is l."2. The 
PSF is expected to be 2". However, I assume that the effective size of the PSF is increased 
to 3" by dithering (see § 4.2). These figures imply that the background is ~ 45 e - s -1 . 

For the great majority of events detected to date, the source star is V > 20. For 
illustration, I consider two V = 20 stars, a main-sequence star (V — L = 0.3) and a 
clump giant (V — L = 2.5). I assume that 1/uJy corresponds to L = 21.1. For the main 
sequence star, I find that the total exposure time required to reach a* = 1% precision is 
T exp = 39(/3 2 + 0.04/5) hr. For the clump giant, the time required is T exp = 0.6(/5 2 + 0.33/3) hr. 
Thus, for main sequence stars with f3 ^ 0.4, the total satellite time is ~ 2T exp ^ 12 hours, 
while for a red giant, the time required is less than one hour. Clearly, the latter is to be 
preferred. To date, unfortunately, only one LMC clump giant event has been published 
(Alcock et al. 1998a). 



4.4. Backgrounds 

As I have previously discussed (Gould 1998b), an asymmetric light curve can be the 
result of a binary lens or a binary source. I made a rough estimate that ~ 20% of events 
could be affected by such backgrounds which would lead to a spurious measurement of 
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Au x . There are no backgrounds that would mimic a shift in the peak time as seen by the 
satellite relative to the Earth, assuming that the peaks were well resolved. Nevertheless, 
the method proposed here is to determine the peak from a pair of symmetrically placed 
observations. Thus, the same asymmetry that produces a spurious detection of 7®, could in 
principle produce a slight shift in to,5- However, 7 is measured from the "wings" of the light 
curve, (t — t )/t e ~ ±1, (Gould 1998b), while the peak is determined from observations 
at (t - t )/t e ~ ±(2/3) 1 / 2 /3. Since dF/d-f oc 0.5(t - t ) 3 while dF/dt cx (t - t )\ the 
effect of any background on the to determination will be ~ (3 2 /3 smaller than on the 7 
determination. That is, it will most likely be negligible. Thus, it would seem advisable to 
push the satellite observations so that o"A Uy < fA^- If the asymmetry detected from the 
ground is truly due to parallax, one should also be able to detect a time delay provided that 
the observations are sensitive enough. Thus, the satellite observations provide a partial 
check on the reality of the ground-based detection of parallax. 

5. Other Lines of Sight 

In this paper I have focused attention mainly on the LMC, partly for simplicity and 
partly because I consider it to be the most interesting line of sight scientifically. However, 
there are two other lines of sight for which one might want to obtain satellite parallaxes: 
the SMC and Galactic bulge. 

The scientific question regarding events detected toward the SMC is similar to that for 
LMC events: are the lenses in the halo or in the Magellanic Clouds? The major difference 
between the SMC and LMC is that the LMC lies almost exactly at the ecliptic pole, 
whereas the SMC lies about 25° from the pole. This difference in turn has two implications. 
First, the equations describing parallax become more complicated. See, for example, Gaudi 
& Gould (1997). The "parallax ellipse" becomes more flattened, making the parallax effects 
smaller and so more difficult to measure. However, since the flattening is only by a factor 
cos 25° ~ 0.9, the effect is quite minor and can be ignored for our purposes. Second, SIRTF 
cannot observe the SMC for the full year as it can the LMC. The telescope can only point 
between 80° and 120° from the Sun. Hence the SMC is only observable for ~ 65% of the 
year. This is not a major problem, but combined with the fact that the SMC is smaller 
than the LMC (and so has fewer source stars to monitor) it does mean that it will provide 
fewer opportunities for parallax measurements. 

The bulge is qualitatively different. First it lies very close to the ecliptic which implies 
that it can be observed only for two 40 day intervals during the year. Second, the parallax 
ellipse is highly flattened. For example, Baade's Window lies only 6° from the ecliptic so 
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the ellipse is flattened by a factor of 10. This flattening will reduce the size of the parallax 
asymmetry by a factor of 10 in the summer and winter, and will reduce the time delay 
between the Earth and satellite by a factor 10 in the spring and fall (Gaudi & Gould 1997). 
Since the bulge is observable from SIRTF only during the spring and fall, it is the latter 
effect that is relevant. Thus, equation fl2lD will be virtually unaffected, but equation (|26| ) 
will be increased by a factor of 10. This degradation of the S/N is partially mitigated by 
the fact that bulge events tend to be shorter (t e ~ 10 days rather than 40 days), but it 
is exacerbated by almost as large a factor because most of the lenses are expected to be 
in the bulge implying that v = {d os /d\ s )v ~ 800 km s -1 (rather than 275 km s -1 ). Thus, 
the photometric precision required is increased by approximately a factor of 7. From the 
standpoint of S/N, this is not a major problem. The bulge is about 6 times closer than the 
LMC, so clump giants are about 36 times brighter. Hence, the exposure time is formally 
only T exp = 0.3/3 hours. However, whether systematics will compromise photometry at the 
required 0.14% level remains an open question. 

From a scientific standpoint, bulge parallaxes can address two principal questions. 
First, where are the lenses? The conventional wisdom is that most are in bulge. However, 
the same conventional wisdom predicts a much lower optical depth and many fewer short 
events than are actually observed. One would like some experimental confirmation of this 
wisdom. Second, what is the mass spectrum of the lenses? Again, the conventional wisdom 
is that the lensing events are due to normal stars in the bulge (and secondarily the disk) 
along the line of sight. However, the observed mass spectrum of bulge stars (Holtzman 
et al. 1998) does not seem to be able to explain the observed distribution of time scales 
(Han & Gould 1996). Han & Gould (1995) showed that parallax measurements could help 
constrain both the location and the mass spectrum of the lenses. 

Even if the time differences At are initially too small to measure, the situation will 
gradually improve with time because the Earth-satellite distance d will gradually grow. 
This tends to increase Ato in two distinct ways. First, of course, as the satellite gets further 
from the Earth, it takes longer for a lens to move from one to the other (see eq. 26). In 



addition, the satellite is constrained to observe the bulge not when it is spring or fall on 
Earth, but when the Sun is near the vernal or autumnal equinox as seen from the satellite. 
As the satellite moves farther from the Earth, the Earth-satellite separation vector becomes 
less closely aligned with the direction of the bulge during these critical times that the bulge 
is observable. Unfortunately, these same changes also make it more difficult to measure 
the parallax asymmetry from the ground. Nevertheless, if it initially proves too difficult to 
measure bulge parallaxes, the situation should be reviewed periodically in light of ongoing 
experience. 
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6. Conclusions 

The method previously developed for measuring microlens parallaxes by directly 
comparing Earth-based and satellite-based photometry will not work for SIRTF. The old 
method requires that both sets of observations be done in the same band in order to remove 
degeneracies in the parallax solution. The SIRTF L band detector response cannot be 
mimicked from the ground because of atmospheric absorption and high background. 

However, by combining ground-based measurements of the "parallax asymmetry" 7 of 
a lensing event (due to the Earth's acceleration) with the observed difference At in the 
peak times of the event as seen from the Earth and SIRTF, it is possible to measure the 
parallaxes of microlensing events seen toward the LMC. The parallax yields the projected 
velocity of the lens, v = (d os /d\ s )v and so would reveal whether the lenses are in the 
Galactic halo or in the LMC itself. 

The ground-based measurement requires ~ 1% photometry about once per day 
for about 2.5 Einstein crossing times, i.e., ~ 100 days. The space-based measurement 
requires two observations each with about 1% precision, plus one additional lower-quality 
measurement. The two 1% measurements should be spaced approximately symmetrically 
about the peak of the event, while the lower-quality measurement is needed to constrain 
the baseline. For a main-sequence star (V = 20, V — L = 0.3), the total satellite time 
required is ~ 12 hours (f3/0A)(d/0.2 AU)^ 1 , where (3 is the impact parameter and d is the 
Earth-satellite distance. For a clump giant (V = 20, V — L = 2.5), less than 1 hour is 
required. 

Events must be alerted in real time and an improvement in the current magnification 
threshold (A ~ 1.6) for the alerts would be helpful but is not critical. Photometry should be 
carried out using image differencing. Backgrounds due to binary lenses and binary sources 
are minor but not completely negligible. 

Parallax measurements are also possible for SMC events. The major difference from 
the LMC is that the SMC can be observed for only 70% of the year. Parallax measurements 
of Galactic bulge microlensing events are substantially different. Formally, the telescope 
time requirements are less severe than toward the LMC primarily because the sources are 
substantially brighter. However, because the bulge lies near the ecliptic, the size of the 
parallax effect is substantially smaller than toward the LMC, and this may mean that small 
systematic errors will compromise the measurements. 
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Fig. 1. — Normalized errors, R {bold), S (solid), and T (dashed), for the measurements of 
to, j3, and 7. These are respectively the time of maximum, the impact parameter, and the 
asymmetry parameter. Each is plotted as a function of the impact parameter, f3. The actual 
error is given by, for example, cr 7 = Sa /(Nt e / 'day) 1 / 2 where cr Q is the factional flux error 
for an individual measurement, and (Nt e / day) is the number of measurements per Einstein 
crossing time. These curves assume that light-curve measurements begin when the source 
enters the Einstein ring (m = 1.0) and end at t = to + 1.5t e . 



Fig. 2. — Ratio / of the actual error in Au x (taking full account of all the covariances 
among the 13 parameters of the full fit) to the naive error given by equation (|23"D (together 
with the approximations that the error in Au x is equal to the error in At /t e , and that 
the latter is dominated by the error in to t s/t e ). This ratio is shown for various values of 
the impact parameter f3 as a function of St. Here, St is the time difference between the 
peak of the event and the midpoint of the two observations by the satellite. If the two 
observations are symmetric about the peak (St = 0), then 1 < / < 1.36. However, the 
errors can increase dramatically for non-symmetric observations, particularly for small (3. 
Symmetric observations yield more precise results because the uncertainties in t e , (3, and 
F — B (which are all relatively large) are then decoupled from the uncertainties in t . The 
curves shown in the figure all assume an Earth-satellite separation d = 0.2 AU and Einstein 
crossing time t e = 40 days. However, for other values of t e <; yr/27r and d <C AU, the curves 
are qualitatively similar. 



Fig. 3. — Normalized error S for the asymmetry parameter 7 which is used to determine 
Aiiy. See equations (|T2|), (|I~9|), and (pT|). The lower curve assumes that intensive follow-up 
observations begin when the enters the Einstein ring (ui n i t = 1) and is the same as solid 
curve in Figure 0. This curve is assumed in all calculations in the paper. The upper curve 
(i*imt = 0.75) reflects the present capability of the MACHO collaboration alert program. If 
this capability cannot be improved, the errors would increase substantially for (3 ^ 0.45. 



